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An external magnetic field, H , applied parallel to a quasi two-dimensional system modifies quan- 
titatively and qualitatively the density of states. Using a self-consistent numerical approach, we 
study how this affects the entropy, S, the free energy, F, and the magnetization, M, for different 
sheet carrier concentrations, N s . As a prototype system we employ III-V double quantum wells. 
We find that although M is mainly in the opposite direction of H, the system is not linear. Sur- 
prisingly dM/dH swings between negative and positive values, i.e., we predict an entirely orbital 
diamagnetic to paramagnetic fluctuation. This phenomenon is important compared to the ideal de 
Haas-van Alphen effect i.e. the corresponding phenomenon under perpendicular magnetic field. 

PACS numbers: 75.20.-g, 75.75.+a, 75.60.Ej 



The scrutiny of quasi two-dimensional (2D) carriers 
under magnetic field has a long, fruitful history. Empha- 
sis was given to the configuration where the magnetic 
field, H, is applied perpendicularly to the quasi 2D sys- 
tem, leading to quantization of the free in-plane motion 
into Landau levels In this configuration the integer 
quantum Hall effect [2j was discovered and the fractional 
quantum Hall effect was observed Q and explained in 
terms of quasiparticles with fractional charge [4]. Be- 
sides, varying H, oscillations of the magnetic susceptibil- 
ity (the de Haas - van Alphen effect [5() as well as oscil- 
lations of the conductivity (the Shubnikov-de Haas effect 
Q) were observed. The heart of the system, the density 
of states (DOS) can be directly probed by measurements 
of a thermodynamic quantity like the magnetization, M. 
Experimental evidence of the ideal de Haas-van Alphen 
effect in a quasi 2D carrier system was only found recently 
by Wilde et al. [7|] who measured the magnetization oscil- 
lations of high mobility 2D electrons in modulation-doped 
AlGaAs/GaAs heterostructures. In the sample with the 
highest oscillation amplitude they observed discontinu- 
ous jumps in M with peak-to-peak amplitude of two ef- 
fective Bohr magnetons (fi* B = <^|qr, m* is the effective 
mass) per electron, in agreement with old Peierls pre- 
diction [§] . Numerical simulations 0] assuming no states 
between the Landau levels (no "background" DOS) could 
model these jumps quantitatively, but measurements on 
samples with lower mobility revealed a finite background 
DOS. Indeed in earlier studies [9( a considerable DOS 
between the broadened Landau levels was necessary to 
simulate the experimental results. 

This article is devoted to the entropy, S, the free en- 
ergy, F, the magnetization, M, and the magnetic sus- 
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ceptibility, \ m — dM/dH, of quasi 2D carriers under 
in-plane magnetic field, employing a self-consistent en- 
velope function approach. In particular, we show that 
although M remains basically in the opposite direction 
of H, the system is highly non-linear. As a result, the 
magnetic susceptibility oscillates between negative and 
positive values. This is the first prediction of a purely 
orbital diamagnetic to paramagnetic fluctuation. The ef- 
fect, ignored by the community up to now, is important 
compared to the ideal de Haas- van Alphen effect, the cor- 
responding phenomenon under perpendicular magnetic 
field. We hope that this article will also help the inter- 
pretation of magnetization measurements under tilted H, 
i.e., whenever an in-plane component of H exists. 

When a quasi 2D system is subjected to an in-plane 
-or even tilted- magnetic field, the charming concept of 
Landau levels must be revised, because carriers move un- 
der the competing influence of the Lorentz force and the 
force due to the quantum well (QW) confining potential. 
The equal-energy surfaces (To| or equivalently the den- 
sity of states [ill EH are qualitatively and quantitatively 
modified because the spatial and the magnetic confine- 
ment compete. Generally, a proper treatment involves 
self-consistent computation [l(J, 13, IH of the energy dis- 
persion, Ei ya (k x ), where i is the subband index, a denotes 
the spin and k x is the in-plane wave vector perpendicu- 
lar to the external in-plane magnetic field (applied along 
y), H. The envelope functions along the "growth" z- 
axis depend on k x i.e., ipi ltr ,k*,k v (r) * Ci,cr,k x {z)e %H « x e %H * v . 
The consequences of this modification were initially real- 
ized in transport [3] experiments. This amendment also 
influences the character of plasmons in single (l5j and 
double 11611 Q Ws. The TV-type kink was theoretically 
predicted [17| an d recently verified in photoluminescence 
experiments [18]. The impact of the DOS modification 
on various properties of dilute-magnetic-semiconductor 
(DMS) single QWs was studied lately pl-lilf. Thus, it 
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seems that the parallel or titled configuration offers new 
avenues to explore. A compact formula for the DOS of a 
quasi 2D system under in-plane H exists (l 31 ] : 
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It is implied that the QWs are along the z-axis and H 
is applied along the y-axis. 6 is the step function, A 
is the xy-axea of the structure. Ei^ a {k x ) are the spin- 
dependent xz-plane eigenenergies which gen erally must 
be self-consistently calculated [13, EEIlE El- Equa- 
tion |T]) is valid for any kind of interplay between spatial 
and magnetic confinement. The k x dependence in Eq. ([T]) 
increases the numerical cost by a factor of 10 2 — 10 3 
in many cases; hence it is sometimes overlooked, al- 
though this is only justified for narrow single QWs or 
for H — > 0. With the existing computational power, 
such a compromise is not needed. In the limit H — > 0, 
Eq. |T]) converges to the staircase shape with the fa- 
mous step for each spin. The opposite asymp- 
totic limit of Eq. fl} is that of a simple saddle point, 
where the DOS diverges logarithmically [ll|. The DOS 
modification drasticall y a ffects the physical properties 
13, EL E3, E3, EI EE 13 E3, E3, El; models which ig- 
nore it can only be applied to narrow single QWs or for 
H — ► 0. For completeness, we note that in Eq. dJ) disor- 
der is ignored. With the progress of the epitaxial tech- 
niques it is fairly small in well-prepared III-V structures. 
Disorder will induce some broadening of the subbands. 

The total population, N, the internal energy, U, the 
entropy [2(J, S, and the free energy, F, are given by: 



Ho is the magnetic permeability of free space. To calcu- 
late M we have to keep the temperature, T, as well as 
N constant e.g. assuming that all dopants arc ionized. 
Here T — 4.2 K. Hopefully, below we list all the remaining 
symbols used in the present article: Ni are the "sheet" 
(1) (measured e.g. in cm -2 ) subband concentrations, while 
N s = N/A is the "sheet" (measured e.g. in cm" 2 ) to- 
tal concentration. We use Ei{k x ) for the subband en- 
ergy dispersions i.e. in the present article we ignore spin- 
splitting which has been treated in detail elsewhere 19]. 
Finally, for H = 0, the "the symmetric -asymmetric gap" , 
AS AS = E x {k x = 0) - E (k x = 0). 
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F = U -TS. 
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fo(£) is the Fermi-Dirac distribution function. Since 
p(£ ) oc A, it follows that N, U, S as well as F are pro- 
portional to A. However, the magnetization, 
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where V is the structure's volume, is independent of A. 
To have the usual units in Tesla, we symbolize B = HoH, 
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FIG. 1: (Color online) QW potential energy profiles, for H = 
0. (a') N s = 1.85 x 10 11 cm" 2 , ASAS = 4.24 meV. (/?') 
iV s = 3.70 x 10 11 cm" 2 , ASAS = 4.01 meV. (7') iV s = 5.55 
x 10 11 cm' 2 , ASAS = 3.79 meV. 

As a prototype system, we choose GaAs/(Al,Ga)As 
double QWs, a bilayer system with well defined 
"symmetric-asymmetric gap" for H — and well-known 
material parameters. Magnetization measurements un- 
der perpendicular magnetic field of such a bilayer 2D 
system with strong coupling between the two QWs can 
be found elsewhere [2lj. To facilitate the reader we pro- 
vide in Fig. Q] the self-consistent potential energy pro- 
files, for H = 0, of the various double QWs employed in 
the present article. Two (left and right) 50 nm spacers 
separate the S— doped layers from the double QW. The 
total double QW width is 22.7 nm including the inter- 
nal barrier of 2.5 nm. For simplicity we take A = 1 m 2 . 
Augmenting the <5-doping, we vary N i: N s and ASAS, 
distinguishing three cases: (a') Nq — 1.51 x 10 11 cm" 2 
and Ni = 0.34 x 10 11 cm" 2 , i.e. N s = 1.85 x 10 11 
cm- 2 , while ASAS = 4.24 meV. (/?') N = 2.41 x 10 11 



cm" 2 and N r = 1.29 



x 



i.e. N. = 3.70 x 



10 11 cm" 2 , while ASAS = 4.01 meV. (7') N„ = 3.30 x 



2 

" 2 and Ni = 2.25 x 10 11 cm" 2 , i.e. N s = 5.55 x 
2 , while ASAS = 3.79 meV. Fig. [2] depicts N , 



10" cm 
10 11 cm 

iVi and N s as functions of figH for all cases. We notice 
that N s (a') : N s (/3') : N s (j') = 1 : 2 : 3. The depop- 
ulation of Ei induced by the DOS modification occurs 
approximately at 6 T, 10 T, and 12.5 T, respectively. 
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FIG. 2: (Color online) The sheet subband concentrations No, 
Ni and the sheet total concentration N s as functions of the ex- 
ternal in-plane magnetic field, fioH , for cases a' (solid lines), 
13' (dashed lines), 7' (dash-dotted lines). 

Fig. [3] depicts U, —TS and F, as functions of HoH, for 
all cases. At T — 4.2 K, \F\ w |U] > | - T5|. Since 
AT is kept constant in each case, we expect that \U\ will 
decrease whenever H induces "flattening" of the occupied 
subbands i.e., expansion of the occupied parts to higher 
\k x \, since this leads to occupied energies with smaller 
\£\. The gradual increase of \F\ from (a') to ((3') and 
(7') mirrors the increase of the population. To facilitate 
the reader we provide in Figs. [?] zooms of the energy 
dispersion of case (a 1 ) for characteristic values of /iqH. 
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FIG. 4: (Color online) Zooms of the energy dispersion, Ei(k x ) 
(i — 0, 1), for characteristic values of ^lqH, for case (a'). Here 
the Fermi energy is identified with zero. 
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FIG. 3: (Color online) The internal energy, U, the product 
—TS and the free energy, F, as functions of the external in- 
plane magnetic field, fioH, for all cases (a , j3 , 7'). On this 
scale, F fa U, —TS is negligible. 

Fig. [5] depicts the entropy S as a function of fj,oH for 
all cases (a', /?', 7'). According to Eqf5]and Eq. [3J since 
for each case N is constant, S 1 is sensitive to the changes 
of ln[fo(£)]- At T = 4.2 K, these changes only occur in 
a short region around the Fermi energy, Ep. In other 
words, S reads the modification of the energy dispersion 
around Ef = 0. For case (a 1 ), for (iqH = 0, the bottom 
of Ei(k x ) is relatively close to Ep = 0. From to 6 T, 
S is continuously falling due to the continuous depopu- 



lation process of E\{k x ). The minimum of S occurs at 
6 T because there the gradual DOS modification com- 
pletely depopulates the 1st excited subband, cf. Fig. [2] 
and Fig. 2J In crude language this means less anarchy, 
less entropy. From 6 T to 10.75 T, S is continuously 
increasing because around Ep = 0, Eo(k x ) in the range 
\k x \ ~ moves upward and finally the populated Eo(k x ) 
is divided into two parts, cf. Fig. [4] The maximum of S 
occurs at 10.75 T because there Eo(k x ) splits completely 
into two occupied parts, around k x — the ground state 
subband ceases to be occupied. Thus, the cohesion of 
the occupied Eo(k x ) is lost around 10.75 T. Hence, more 
anarchy, more entropy. From 11 T to 20 T, S does not 
change much, because the main effect of increasing i^qH 
is to move the two Eo(k x ) minima continuously apart, 
cf. Fig. [U The behavior of S, in cases (/?', 7') can be 
readily explained in analogous manner. As one could 
maybe imagine, increasing the system's magnitude, the 
minimum of entropy, S m i n , is increasing: at 6 T with 
S mm = 2.2 xlO 10 eV/K, at 10 T with S min = 2.4 xlO 10 
eV/K, at 12.5 T with S mi „ = 2.5 xlO 10 eV/K. 

Fig. [6] shows the magnetization M as a function of 
HoH for all cases (a', /?', 7'). We observe that the DOS 
modification induces an oscillation of M; it is between 
« -3 A/m and « 0.5 A/m for case (a'), between -7 
A/m and « 2 A/m for case (/?'), and between w -9 A/m 



4 



— P' 






— y' 



















U ^ 1 — i — , — i — , — i — , — i — , — i — , — i — , — i — , — i — , — i , 

2 4 6 8 10 12 14 16 18 20 

f H(J) 

FIG. 5: (Color online) The entropy, S, as a function of noH, 
for each one of the cases (a' , /3' , 7'). 

and pa 2 A/m for case (7'). The reader may observe 
that the magnetic susceptibility, Xm — dM/dH, swings 
between negative and positive values, thus Fig. [5] shows a 
totally orbital diamagnetic to paramagnetic fluctuation. 
For example, for the case (7') the fluctuation of M of the 
order of 10 A/m, is translated to approximately ^ of the 
ideal de Haas-van Alphen effect i.e. the magnetization 
step of two effective Bohr magnetons per electron in the 
perpendicular configuration 0, Q ■ 
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FIG. 6: (Color online) The magnetization M) as a function 
of fJ.oH, for each one of the cases («', /3', 7'). 

In summary, the main message of this article is that 
an in-plane magnetic field causes a considerable effect 
in the magnetization of quasi two-dimensional carriers 
which has been ignored by the community up to now. 
The magnetic susceptibility swings from negative to pos- 
itive values i.e. we have proved that an entirely orbital 
diamagnetic to paramagnetic fluctuation exists. We con- 
jecture that the in-plane component of a tilted magnetic 
field will bring about similar effects, hence care must be 
taken for the interpretation of such magnetization mea- 
surements. The entropy of quasi 2D carriers at low T 
was also calculated and interpreted. 

I thank Dr. I. M. A. Bominaar-Silkens and Dr. U. 
Zeitler for useful discussions. 
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